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We study the A'^-dependence of the thermodynamical variables and the dynamical behavior of the 
well-known Hamiltonian Mean Field model. Microcanonical analysis revealed a thermodynamic 
limit which defers from the a priory traditional assumption of the A'^-dependence of the coupling 
constant g as ^ 1/N: to tend N ^ oo keeping constant E/N^ and g/N, prescription which guar- 
antees the extensivity of the Boltzmann entropy. The analysis of dynamics leads to approximate 
the time evolution of the magnetization density m by means of a Langevin equation with multi- 
plicative noise. This equation leads to a Fokker-Planck's equation which is Af-independent when 
the time variable is scaled by the A'-dependent time constant Tmac = \/ IN/g, which represents the 
characteristic time scale for the dynamical evolution of the macroscopic observables derived from 
the magnetization density. This results explains the origin of the slow relaxation regimen observed 
in microcanonical numerical computations of dynamics of this model system. Connection with the 
system self-similarity is suggested. 
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I. INTRODUCTION 

The present effort is devoted to study the A^- 
dependence in the thermodynamical variables and dy- 
namical behavior of the well-known the Hamiltonian 
mean field model [111110,110. This model describes 
a system of A^ planar classical spins interacting through 
an infinite-range potential : 

Hn {0, L; I, g) = Kn {L; I) + Vn {0; g) 



1 ^ 



(1) 



where 6i is the i-th angle and Li the conjugate variable 
representing the angular momentum, / is the moment of 
inertia of the rotator and g the coupling constant. This 
is an inertial version of the ferromagnetic X-Y model , 
which interaction is not restricted to first neighbors but 
is extended to all couples of spins. 

Traditionally the coupling constant g is suppose to be 
A^-dependent: 
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assumption which makes LL formally extensive {V (x N 
when A^ ^ oo), since the energy remains non-additive 
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and the system can not be trivially divided in inde- 
pendent subsystems. A fundamental observable of this 
model system is the magnetization density m, which is 
given by: 



1 ^ 



A^ 



(3) 



1=1 



where rrife — [cos^fe, sin^^] is the intrinsic magnetization 
of the fc-th particle. 

The canonical analysis of the model predicts a second- 
order phase transition from a low-energy ferromagnetic 
phase with magnetization m « 1 to a high-energy para- 
magnetic phase with m « 0. The dependence of the en- 
ergy density U — E/N on the temperature T, the caloric 
curve, is given by: 



{T)] 



(4) 



{U and T in units of 7) and this is shown at the Fig.Q). 
The critical point is at the energy Uc = 0.75 corresponds 
to the critical temperature of Tc = 0.5 0. 

The dynamical behavior of the HMF has been exten- 
sively investigated. Microcanonical simulations of dy- 
namics are in general in good agreement with the canon- 
ical solution, except for a region below Uc, where some 
anomalous phenomena have been observe, such as en- 
semble inequivalence and the slow relaxation regimen to- 
wards the Boltzmann-Gibbs' equilibrium il0,H|^. 

This ensemble inequivalence, characterized by the exis- 
tence of a negative specific heat, appears as a dynamical 
feature. The origin of this anomalie is explained by the 
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existence of certain quasi-stationary state |^|^ which is 
not the canonical one and depends on the initial con- 
ditions in the microcanonical numerical computation of 
the dynamics. In this case, the system needs a long time 
to relax to the canonical equilibrium state because of 
the duration of this metastable state has a linear depen- 
dence with N. This fact has been interpreted by Tsallis 
as the non-commutativity of the thermodynamic limit 
(N — !■ oo) with the infinite time limit ( t — > oo): when 
the first is performed before the second, the system will 
not relax to the Boltzmann-Gibbs' equilibrium. 

Recent studies have shown that the slow-relaxation 
regime at long times is clearly revealed by numerical re- 
alizations of the model, but no traces of quasistationarity 
are found during the earlier stages of the evolution. This 
fact suggested the nonergodic properties of this system 
in the short-time range, which could make the standard 
statistical description unsuitable ^6]. 

However, it there could be a mistake in the under- 
standing of the macroscopic description of this model. 
There are many studies in long-range Hamiltonian sys- 
tem where the validity of the thermodynamical descrip- 
tion in the thermodynamic limit is intimately related to 
the scaling with N of the thermodynamic variables and 
potentials. The problems with the nonextensive nature 
of those systems are avoided using the Kac prescription 
[gj in which the coupling constants are scaled by some 
power of N in order to deal with an extensive total en- 
ergy E. It is well-known that this kind of assumption 
affected the time scale of the system evolution, so that, 
its application should be taken with care. 



II. MICROCANONICAL ANALYSIS 



W = We are interested in describing the large N 

limit, where W is given by: 



^"(H)' /rf'mK + 2C/-l]^^/(m;Ar), 

(7) 

where it was introduced the dimensionless parameter U: 



and the fimction / (m; N): 
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fc=i 
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We estimate the function / (m; N) for large N by using 
the steepest decent method. Expressing the delta func- 
tion by mean of the Fourier's integral representation, the 
Eq.© is rewritten as follows: 



2 exp [iNk ■ m] 

(27r) 



lo 



N 



where k = k-fix, x e R^, and Io{z) is the modified 
Bessel function of zero order. It is easy to obtain the 
main contribution of the integral when N ^ oo maxi- 
mizing the integral function via the parameter x. That 
is: 



M. Antoni, H. Hinrichsen and S. Ruffo carried out the 
microcanonical analysis of the HMF in the ref. 9]. How- 
ever, in their work they considered the equivalence of the 
microcanonical with the canonical ensemble in the ther- 
modynamic limit, which depends on the validity of the 
canonical solution. In this section we will perform this 
analysis working directly on the microcanonical ensem- 
ble. 

The microcanonical states density is calculated as 
follows: 



1 f d'^ f)d" T 
n{E,N;I,g)^— / -^^5 [E - {9, L;I, g)] 
■ J [2nn) 



The integration by d L yields: 



(5) 



/ (m; N) ~ exp {-N [xm - In Iq (x)] 



■In 



/ NV m{x) d , ' 

t:— i-'fn (x) 

\2tt J X dx ^ ' 



where x = |x| is related with m = |m| by: 



m = m (x) 



hjx) 
lo (x) ■ 



(10) 



(11) 



It is necessary to say that the validity of the steepest 
decent procedure is ensured by the positivity of the ar- 
gument of the logarithmic term in the Ea. (|10|) . 



1 /27r/\^" f d^e , ,„ ^nlTV—i For continuous variables, W is only well-defined after a coarsed 

NW f ^ ") \ fi2 ) / To Tn ( ' 3)\ ^ ■ grained partition of phase space, which is the reason why it is 

\2/^' v""/ considered an small energy constant <5e in order to make W di- 

l^j mensionless. However, during the thermodynamic limit —* oo, 

The microcanonical accessible volume W is expressed as the selection of 5e does not matter if it is small. 
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The main contribution to the entropy per particle 
s {U,N;I,g) — InW/N for N large can be obtained by 
using again the steepest decent method as follows: 



+max 2 t™'^ + — xm (x) + In Iq (x) ^ , (12) 
where k — 2U — 1. The maximization yields: 



m (x) 



rr? (x) — K 



w? (x) + K dx 



m{x)=0, (13) 



which represents the energy dependence of magnetization 
density for the bulk matter: 



TO — m{x) 
U=^m{x) /x+^[l-m^ (x)] 

and TO = if U > Ur ^ 0.75 



with X G [0, oo) 



(14) 

The caloric curve (0)) is obtained from the canonical pa- 
rameter /3: 



1 _ dso{U,N: 
^ T dU 



1,9) 



x/m{x) with 0<U < 0.75 
1/(2C/-1) if [/>0.75 



(15) 



This is the same result found by M. Antoni, H. Hinrich- 
sen and S. Ruffo in the ref.'^. Therefore, the subsequent 
analysis will be equivalent to their work. The micro- 
canonical analysis confirms the existence of a second- 
order phase transition at the critical energy Uc = 0.75 
with a critical temperature of Tc = 0.5. 

It is interesting to note that although the main term 
of the entropy per particle is A/^-independent, there is an 
additional contribution which is iV-dependent and makes 
the Boltzmann entropy nonextensive for an arbitrary A'^- 
dependence of the coupling constant g. Although the se- 
lection of the A^-behavior of the coupling constant does 
not change the thermodynamical picture in the micro- 
canonical ensemble, it affects the extensive or nonexten- 
sive character of the energy in the thermodynamic limit, 
as well as the temporal scale in the dynamical behavior 
of this model system (see in the next section). 

In order to deal with an extensive Boltzmann entropy 
is necessary to assume the following thermodynamic limit 
for the Hamiltonian mean field model: 



E g 

N oo, keeping constant and — . (16) 

This assunition defers from the one used by other investi- 
gators, where the A^-dependence of the coupling constant 
chosen in order to deal with an extensive energy. 



FIG. 1: Microcanonical description of the HMF model in ther- 
modynamic limit. Here it is clearly shown a second-order 
phase transition at (7c = 0.75, where Tc = 0.5 and m — 0. 



III. THE ORIGIN OF THE SLOW 
RELAXATION REGIMEN 

According to the Hamiltonian of the model, its motion 
equations are given by: 



h — jLk, 



Lk = gN {niy cos Ok - nix sin 6k) 



(17) 



where fc = 1, 2, ...A^ and it was used the expression given 
in Eq.Q. It is convenient to introduce the following 
characteristic units for time and momentum: 



To 



— and Lq = ^/igN. 
gN 



(18) 



Note that, according to the thermodynamic limit H16|l . 
the characteristic time tq decreases with the A^ increas- 
ing as To (X 1/A", while the characteristic momentum Lq 
increases as Lq oc N . However, when it is assumed the 
A-dependence for the coupling constant, the char- 
acteristic time To is A-independent. The time scale of 
To gives the characteristic time scale for the evolution 
of each rotator, so that, this time units is relevant to 
microscopic level. In the numerical computation of the 
dynamics this characteristic time scale should be taken 
into account in order to ensure a great accuracy of the 
energy conservation. 

Using the above characteristic units, as well as intro- 
ducing the tridimensional vectors rrife = (cos Ok , sin 0fe , 0) 
and Ofe = (0,0, Lfe), the motion equations (fTTjl are 
rephrased as follows: 



mfc = Ofe X niA;, Ofc = nife X m. 



(19) 
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The magnetization density obeys the following dynamical 
equation: 



m=T = Jlxm + Am, 



(20) 



where T is the velocity of change of the magnetization 
density, fJ, the angular velocity which characterizes the 
dynamical behavior of the m orientation, while A char- 
acterizes the behavior of its modulus. These quantities 
are given by: 



1 ^ 



(21) 



fe=i 



1 ^ 1 ^ 

k=l k=l 

Since x rill's do not have a definite sign, the magne- 
tization velocity of change Y is a fast fluctuating quan- 
tity with a vanishing short time average (a time average 
along temporal interval comparable with the characteris- 
tic time of rotators evolution tq) and standard deviation 
ctt '^1/ VN. This fact provokes that the short time av- 
erage of the magnetization m varies slowly in comparison 
with the characteristic time scale of rotators evolution. 
Thus, there is two time scales in the dynamical evolution 
of the system: the microscopic time scale, Tmic = tq, 
which is given by the characteristic time of the rotators 
evolution, and the macroscopic time scale, Tmac, which is 
given by the characteristic time of the short time average 
of magnetization density evolution, being this last con- 
siderably greater than the first time scale, Tmac 3> Tmic ■ 

We are interested in studying the dynamical evolution 
of the macroscopic observables for temporal scale compa- 
rable with Tmac- Since Tmac ^ Tmic, at first approxima- 
tion the dynamical evolution of m could be modeled by 
a Langevin equation with multiplicative noise |lCl : 



m = T ==(Ti (m; U) Ci (t) ei + <T2 {m; U) ^2 {t) 63. (23) 

whose parameters cti (m; U) and 02 (to; U) could be esti- 
mated by using a microcanonical average with fixed mag- 
netization density. This last estimation is a reasonable 
approximation since the magnetization density only ex- 
hibits very small fluctuations around its short time av- 
erage for N large along time intervals comparable with 
Tmic- Here, ei and e2 are unitary vectors in the paral- 
lel and perpendicular directions of m respectively, being 
63 = ei X 62 = (0, 0, 1), the unitary vector along the third 
axis, [t) and ^2 (*) are Gaussian processes: 



(e. {t)tj{t'))^TA,5{t^t') 



(24) 



being the correlation times. The presence of 5ij is 
due to the vanishing of the microcanonical average of 



(Y1Y2) = (Y • eiY • 62) with fixed magnetization den- 
sity, which is easily proved by mean of a procedure ana- 
logue to the one used in the previous section for the deter- 
mination of Boltzmann entropy. As already mentioned, 
the standard deviations of the components of the vector 
Y depend on N as Cj (to; U) ~ l/y/N. 

It is not difficult to understand that the correlation 
times Ti are comparable with the microscopic time scale 
Tmic- the correlation times possess the same temporal 
scale of the characteristic time where the differenciable 
character of Y is perceived, which is precisely character- 
istic time scale of rotators evolution. 

According to the general theory, to this equation is 
associated the Fokker-Planck equation of the form: 



-F-d 
dt ~ 



juje"d^ 



o,e]F 



(25) 



which describes the dynamics of the magnetization den- 
sity distribution function F = F (to, t; U). In this expres- 
sion, the e" 's represent the components of ej . Taking into 
account the A^-dependence of aj, the right hand of the 
equation H25|l decreases as during the N increasing. 
Scaling the time scale as t — > t = Nt', the Fokker-Planck 
equation becomes A^-independent. This means that the 
characteristic time scale for the dynamical evolution of 
F (to, t; U), as well as all those macroscopic observables 
derived from m, have a linear dependence on N during 
the increasing of system size (in tq units). 

This is precisely the result which was observed in the 
numerical computation of the dynamics, where the re- 
laxation time for the dynamical temperature To (<) = 
2K (t) /gN"^ = 2C/ - 1 -t- TO^ (t) grows proportional to A^ 
with the A^ increasing. As already mentioned, the time 
scale of such dynamical behaviors is precisely the macro- 
scopic time scale Tmac- 



ToN 




(26) 



Thus, when the thermodynamil limit (|16|l is taken into 
account, the relaxation time is finite when A^ is tended 
to infinity. Therefore, no anomalous dynamical behav- 
ior will be observed: the thermodynamic limit (A^ — > 00) 
necesary for the ensemble equivalence between the mi- 
crocanonical ensemble with the canonical one, will com- 
mute with the infinite time limit (t 00) necesary for 
the equilibration of the temporal average of the physical 
observables. Thus, the slow relaxation regimen observed 
in the microcanonical numerical computation of dynam- 
ics appears as consequence of the consideration of the 
A^-dependence of the coupling constant The above 
analysis shows that the A^-dependence of the coupling 
constant which ensures the extensivity of the Boltzmann 
entropy is more appropriate than the one which guaran- 
tees the extensivity of the energy. The first allows us to 
carry out a well-defined macroscopic description for the 
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Hamiltonian mean field model, instead, the second leads 
to the existence of dynamical anomalies. 



IV. CONCLUDING REMARKS 

As already showed, the microcanonical analysis re- 
vealed a thermodynamic limit H16(l which defers from 
the traditional assumption of the iV-dependence of the 
coupling constant jSJ- This thermodynamic limit is in- 
timately related with the existence of the following self- 
similarity scaling behavior of the thermodynamical vari- 
ables of the Hamiltonian mean field model: 

N N{a) = aN^ ^ 

E^E{a)=a^E, \ ^ ^ {a) ^ T [W {1) M (27) 
g ^ g{a) = ag, [ ^ v y > j 

I^I{a)^I, J 

where the functional [W, a] defines an exponential self- 
similarity scaling laws for the microcanonical volume W: 

T[W,a]=exp[a\nW], (28) 

being 



W{a) = W [E (a) , N (a) ; g (a) , I («)] , (29) 
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